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We study the stochastic dynamics of coupled states with transition probabilities depending on 
local persistence, this is, the time since a state has changed. When the population has a preference to 
adopt older states the system orders quickly due to the dominance of the old state. When preference 
for new states prevails, the system can show coexistence of states or synchronized collective behavior 
resulting in long ordering times. In this case, the magnetization m{t) of the system oscillates around 
m{t) — 0. Implications for social systems are discussed. 
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Models of two states are commonly used in physics as 
a tool to study the emergence of collective behavior in 
systems from spin interaction to opinion dynamics P- 
13]. In the adoption of traits [3HZ] different aspects have 
been studied including the relevance of the interaction 
topology [HHlQl, social influence [HI |T2| , and mass media 
[I3HI11- When accounting for opinion dynamics, the ma¬ 
jority of models are based on decision rules that consider 
a fraction of the surrounding states, e.g., voter model 
m, threshold model HU, majority rule m, or Sznajd 
model [19]. The timing of the interactions can also af¬ 
fect the behavior of the system at least by two ways: 
the precise sequence of interactions and by the aging of 
states. For example, in epidemic spreading and diffusion, 
the temporal sequence of interactions can slow down the 
spreading process [201123] : in ordering dynamics, state- 
dependent updates can have a qualitative impact on the 
mean time to order [MI129| . Aging in physical systems 
refers to the persistence time, that is, the time spent in a 
given state, and affects the response of the system to an 
external field or perturbation [301 El]. In social systems, 
when individuals make choices they usually relay on their 
own past experience or memory [32H34j . While conser¬ 
vative groups tend to hold ideas in an unaltered form for 
a long time, progressive individuals embrace new opin¬ 
ions, ideas, or a technology and disseminate them with 
more enthusiasm |351 136] . In the competition between 
new and old information, although new information is 
more valuable for exploring and spatial searching 133: 
adopting older strategies can promote cooperation and 
group success [35]. Also in a biological context, aging in 
speciation events has been proposed as a mechanism to 
explain the shape of evolutionary trees [39] . 

Here we analyze how the tendency of agents towards 
the adoption of established vs. novel traits influences 
the macroscopic dynamic and the ordering process. We 
tackle this problem by considering a model in which the 


adoption of states depends on the time span the agents 
have held their current states. 

The model is defined as follows: each agent has a state 
I that can be up (f) or down (/) with age young {y) or old 
(z). Agents can be then in four states , y^, z^, and z-^. 
Young agents become old at a rate that we set to r = 1. 
Then, there are reactions of randomly paired agents: i) 
in young i and old j of opposite opinions, i adopts the 
opinion of j with probability w = 5 + e, otherwise (with 
probability 1 — w), j adopts the opinion of i; ii) in pairs 
of agents with the same age and different opinion, each 
agent has probability | of convincing the other; iii) for 
pairs of agents with the same opinion, nothing happens. 
When an agent adopts an opinion, it goes to the young 
age of the adopted opinion. Neglecting correlations, the 
expectation values of state concentrations evolve accord¬ 
ing to 

y^ = —y^ + (1 — 2w)y^z'^ -j- wy^z^ -h iz^z"*" , 
y'^ = —y^ -j- wy^z'^ + (1 ~ 2w)y^z^ -j- ^z^z"*" , 

( 1 ) 

it = yt _ (1 _ _ l^t^t ^ 

it = yt _ (X — rcjyt^t _ , 

with the normalization yt-|-yt-|-zt-|-zt = 1. Here we use 
ytt and ztt to refer to the fraction of the corresponding 
states occupied by the agents. The parameter e corre¬ 
sponds to the persuasiveness of the agent, e > 0 means 
that agents with older opinions are more persuasive. On 
the contrary, e < 0 corresponds to agents with younger 
opinions been more persuasive. 

The system presents three stationary solutions in the 
relevant range of all four variables being non-negative. 
Two fixed points are the homogeneous solutions Si hav¬ 
ing z"*" = 1 and S 2 having z^ = 1. Here either all opinions 
are down (S'!) or all are up ( 52 ) and old. The homoge- 
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FIG. 1 . Eigenvalues of the stationary solutions of Eq. Q 
as a function of persuasiveness e for: (a) the homogeneous 
solutions Si and S2 having the same eigenvalues, and (b) 
the solution S3. Black solid (gray dashed) lines represent 
real (imaginary) parts of the two complex conjugate eigenval¬ 
ues. Doted-dashed black line represent the third eigenvalue 
(real). The fourth eigenvalue (not shown) for the eigenspace 
in (1,1,1,1) direction is zero due to conservation of normal¬ 
ization. 


neons solutions are stable if e > 0. Non-zero imaginary 
parts of two eigenvalues are obtained for e > 1/4. The 
third fixed point S 3 is an up-down-symmetric solution 
with values 

where A = y^25 -I- 4(e^ -I- 3e). As shown in Figure lllb), 
it is stable if e < 0, thus complementary to the stabil¬ 
ity of the homogeneous solutions. A transition from zero 
to non-zero imaginary parts of two eigenvalues occurs 
when e falls below approximately —0.39. In this regime 
of strongly negative e, the system oscillates when relax¬ 
ing from a perturbation out of the symmetric fixed point 
solution S 3 . This stability scenario is qualitatively main¬ 
tained when r changes. As r —> 0, the point at which 
the non-zero imaginary part of the eigenvalues shows up 
shifts towards e = 0. 

Model with continuous ages. We now move from two 
ages to a continuous age space and introduce an age de¬ 
pendent probability. The model is now described as fol¬ 


lows. Agents can be in one of the two opinions, up or 
down. The state of each agent has associated an age de¬ 
fined as: Ti = t — ti, where t is the current time and ti is 
the time when the current state of agent i was acquired. 
The system evolves by randomly selecting a pair of agents 
that, if they are in different states, with probability Pi^j 
agent agent j adopts the state of agent i, and with prob¬ 
ability 1 —pi^j the contrary happens i copies the state of 
j. If the agents are already in the same state, no change 
is applied. After N updates, time t is increased to t -|- 1. 
We define the probability Pi^j to be dependent on the 
age of the states of both interacting agents as 


1 


( 2 ) 


Initially, each agent has randomly assigned one of the 
two opinions and the initial ages ti are uniformly dis¬ 
tributed proportionally to the system size N. We con¬ 
sider random mixing where each agent is allowed to in¬ 
teract with any other agent. The case a = 0 corresponds 
to an updating probability of Pi^j = 0.5 which leads to 
voter model dynamics m- Large values of the exponent 
{a —> 00 ), correspond to situations in which the agent 
with the initial oldest state is imposing her opinion to 
the others. In the other extreme case (a —> — 00 ), the 
youngest state is imposed. 

For q : —?► 00 the system ends up in the state of the old¬ 
est opinion while for values of a G (—oo,0] the system 
adopts any of the two opinions with equal probability. 
For a S (0,oo) the probability that the system adopts 
the state of the initial oldest opinion grows with increas¬ 
ing a but it tends to 1/2 when N increases. 

Figure shows the probability distribution function of 
the density of states p as a function of a. For a = 0, 
the density of states is homogeneous corresponding to 
an equiprobable distribution of states. For a negative 
but close to zero, the dynamic is concentrated around 
p = 0.5, which corresponds to a configuration where 
the agents alternate between any of the two opinions. 
This situation changes gradually to a more homogeneous 
distribution of states as a becomes more negative. For 
a > 0, the states are concentrated around p = 0 and 
p = 1 showing that the system eventually orders in one 
of the two opinions (the presence of the two peaks is due 
to the alternation in being the oldest opinion during the 
initial conditions). 

Figure|^shows the ordering time Si<i{oi), i.e., the time 
that the system needs to reach a final state where all the 
agents have the same opinion, computed as the median 
of the distribution of ordering times from different sim¬ 
ulations and rescaled to the value S'Ar(a = 0). S'Ar(O) 
increases linearly with the N as it does for the voter 
model PHD]. For values a > 0, SM{a) gets smaller than 
S'jv (0) implying that the system orders faster than in the 
voter model case. There is a transition when a crosses 
zero. For values a 0, S]\[{a) increases very fast with N. 
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FIG. 2. Probability density function (codified as a colormap) 
of the dependence of the density of states with a. The den¬ 
sity of states is calculated by computing the normalized cu¬ 
mulative histogram of the fraction of population in each state 
during a fixed simulation period and averaged over 10 ^ real¬ 
izations. The system size was fixed to = 100. 

This is in agreement with the observed dynamics around 
a = — 1 (see Fig.[^. The extreme values of a correspond 
to the cases where, when confronting two states, the old¬ 
est opinion always induces the change (a = -|-oo) or the 
youngest opinion always induces the change (a = —oo). 
The inset of Fig shows the scaling with system size of 
S'jv(a) in these limits. For a = -|-oo the ordering time 
scales as S'a? ^ with the exponent 7 = 1.2. In the 
other limit, for a = — 00 , the ordering time scales as 
S'tv N exp{bN) with b = 0.009. 

In the regime a < 0, what is the behaviour of the sys¬ 
tem during the long ordering times? Figure |^a) shows 
time series of the magnetization of the system. For a 
negative and sufficiently far from zero, the magnetiza¬ 
tion oscillates around zero. Figure |^b) provides further 
analysis by the autocorrelation functions of the magne¬ 
tization time series. The onset of oscillations is observed 
when a passes a value around —0.5 from above. Fig¬ 
ure l^c) shows frequency w and decay constant 7 ex¬ 
tracted from the autocorrelation functions. These values 
do not exhibit significant dependence on system size. The 
decay constant is maximum at the onset of oscillations, 
i.e. where the frequency w becomes non-zero. Both the 
onset of oscillations and the decay behaviour are captured 
by the basic model, cf. Figure Bb). At the transition 
to non-zero imaginary parts (oscillations), the stability 
of the symmetric fixed point solution (S-^ is maximal, 
meaning that perturbations decay fastest. 

To understand further the dynamic around a = 0 
we define a quantity called the convincingness z. Let 
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FIG. 3. Rescaled ordering time 5'jv(q)/*S'jv(0) versus a for 
different system sizes. Open symbols stand for the median of 
the ordering time normalized to the median of the ordering 
time at a = 0. The horizontal black line has been added for 
visualization purposes. Inset: Scaling of the median of S'w (a) 
in the limits a —>■ 00 (solid symbols) and a —>■ —00 (open 
symbols). Dashed lines fit respectively S'jv(+oo) ~ N'^ with 
7 = 1.2 and Sn{—oo) ~ N exp{bN) with b = 0.009. 




FIG. 4. Oscillations and decay of correlations for the finite- 
size model in the regime a < 0. (a) Excerpts from time series 
of the magnetization for two different values of a, system size 
N = 50000. (b) Autocorrelation functions from magnetiza¬ 
tion time series of length T — 10®, system size N = 50000. (c) 
Oscillation frequency 10 (curves with symbols) and decay co¬ 
efficients 7 (no symbols) extracted from time series under dif¬ 
ferent values of a and N. Gurves for different system sizes A 
are almost indistinguishable. Autocorrelation functions A(r), 
T G [ 0 , 100 ], are considered for a least squares fit of the func¬ 
tional form Afit(t) = exp(— 7t) cos{u}t). 
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S^,S~ S [N] be the two sets of agents with equal opin¬ 
ion within each set and different opinions across sets. We 
define the convincingness of vs. S~ as the probability 
z that the interaction of a uniformly random pair of an 
5'+ agent and an S~ agent leads to adoption of the S'^ 
opinion, 


^ = is-^|-‘is-r‘E E ^4^- (3) 

»GS+ jes- * 3 

In case a < 0, there are competing effects governing the 
dynamics of z. When i S S~^ convinces j S S~, i) the 
set S'"*' gains another member who now has the youngest 
opinion increasing z. ii) The set S~ loses a member j 
with Tj typically larger than average, making opinions of 
S~ members younger on average decreasing z. iii) With 
time advancing, all opinions age by the same additive 
rate. This makes ratios between ages smaller, driving z 
towards 1/2. In the case a ^ — 1, the first effect dom¬ 
inates. Thus, an initial advantage in z is amplihed and 
the system orders quickly. For a « — 1, ordering times 
are large due to dominance of the second and third ef¬ 
fects. In order to verify this idea, we numerically record 
data pairs (z(t), z(t + t) — z{t)) with r = 1. Averaging 
over pairs with the same or similar z(t) values, we ob¬ 
tain (z(t-l-r) — z(t)) as the expected restoring force. The 
corresponding standard deviation is the noise strength at 
this z value. The restoring force for z is linear around 
the equilibrium point z{t) = 0.5 while the noise strength 
is mostly independent of z (see Fig. 1 at Supplemental 
Material). This suggests to picture the dynamics around 
q; = 0 as one-dimensional equilibrium in a hyperbolic 
potential under state-independent additive noise. 

Different real systems display dominance such as in 
the adoption of innovations 00 and alternation as in 
opinion formation dynamics |41II43j or economic cycles 
[44] . As an example, Fig. shows the electoral results 
of the governmental elections for United States, United 
Kingdom, and Canada during several decades [ISj. The 
Lomb periodogram [46] of the binary time series reveals 
the existence of alternation between the political par¬ 
ties, by the presence of prominent peaks well above the 
noise level (shuffling of the data), with periods of 20 — 30 
years in agreement with observations |47j . This period of 
time coincides approximately with the length of a gen¬ 
eration. Different mechanisms have been proposed to 
explain political cycles: electorate disappointment |48j . 
voters mood changes [47] |49|, negativity effect [50| or 
policy inertia m- Our study complements those mech¬ 
anisms and contributes to understand how a continuous 
age state dynamics model with competition between pref¬ 
erence for old versus young opinions leads to alternation 
in the leading opinion. 

Summarizing, we have studied the competition of 
states using a basic model that takes into account the 
aging of the state. The stability analysis of the solutions 
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FIG. 5. Election results for United States, United King¬ 
dom, and Canada. Elected parties are represented by their 
official colors, USA: Republicans (red) and Democrats (blue); 
UK: Conservative (blue), Liberals (yellow), and Labour (red); 
Canada: Liberals (red). Unionist (gray), and Conservatives 
(blue). Bottom panel: Lomb periodograms of the binary time 
series for each country. Dashed lines represents the level of 
noise as the result of shuffling the data 250 times. 


reveals the existence of two stable solutions for positive 
values of the persuasiveness (old opinion prevails) that 
compete for consensus. For large negative persuasive¬ 
ness (young opinion prevails), only one solution is stable 
leading to oscillatory transients. We have extended our 
study to a more detailed continuous age model finding 
that, when confronting two opinions, the final config¬ 
uration where only one opinion survives and the time 
needed to reach it is noticeably sensitive to the age of 
the state through the exponent a of the convincing prob¬ 
ability. The continuous age model exhibits likewise the 
oscillations shown by the basic model for large negative 
persuasiveness as well as similar solutions with a. Our 
study provides an alternative mechanism in the under¬ 
standing of the dynamics of consensus formation and the 
observed alternation between states of different systems. 
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